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Abstract 

Using several relativistic mean field models (such as GMl, GM3, NL3, TMl, FSUGold and lU- 
FSU) as well as the quark-meson coupling model, we calculate the particle fractions, the equation 
of state, the maximum mass and radius of a neutron star within relativistic Hartree approximation. 
In determining the couplings of the isoscalar, vector mesons to the octet baryons, we examine the 
extension of SU(6) spin-flavor symmetry to SU(3) flavor symmetry. Furthermore, we consider the 
strange {a* and 0) mesons, and study how they affect the equation of state. We find that the 



'nI" . equation of state in SU(3) symmetry can sustain a neutron star with mass of (1.8 ~ 2.\)Mq even 

if hyperons exist inside the core. In addition, the strange vector {(j)) meson and the variation of 
baryon structure in matter also play important roles in supporting a massive neutron star. 
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I. INTRODUCTION 

Neutron stars, which comprise hadrons and leptons as remnants of supernovae explosions, 
may be beheved to be cosmological laboratories for dense nuclear matter. However, their 
detailed properties, for instance, the mass, radius and particle fractions in the core of a neu- 
tron star, are not fully understood yet, since the pioneering paper by Baade and Zwicky p| 
and the first discovery of a neutron star by Hewish and Okoye [2|]. Because the observed 
mass and/or radius of a neutron star can provide strong constraints on the equation of state 
(EoS) of dense nuclear matter, many theoretical discussions have been focused on the EoS 
to understand the structure of dense matter. 

The typical mass of neutron stars is known to be around IAMq (M©: the solar mass) |3|. 
The most famous, precisely observed pulsar is the binary pulsar PSR B1913+16 (the Hulse- 
Taylor pulsar) with the mass of 1.4398 ± O.OOO2M0 \4, l5|. However, a few neutron stars 
whose masses are much heavier than IAMq have recently been observed. For example, 
Shapiro delay measurements have indicated that the binary millisecond pulsar PSR J1614- 
2230 has the mass of 1.97 ± O.O4M0 [6|. Then, such heavy neutron stars have attracted a 
lot of interest not only in astrophysics but also in nuclear physics, because of the possibility 
of exotic degrees of freedom, such as quarks, gluons and/or some unusual condensations of 
boson-like matter, in the core. 

Recently, many people have often used relativistic mean-field (RMF) models (or relativis- 
tic Hartree models) including hyperons {Y) to calculate the EoS for a neutron star. However, 
it is quite difficult to explain the heavy neutron stars by such EoSs with the meson-baryon 
coupling constants based on SU(6) (quark model) symmetry, because the degrees of freedom 
of hyperons make the EoS very soft, and thus the possible maximum mass of a neutron star 
is considerably reduced P, 18| . 

In Refs. [9I, [lO|, the properties of a neutron star have been studied in detail within 
relativistic Hartree-Fock (RHF) approximation. In those calculations, we have considered 
not only the tensor couplings of vector mesons to the octet baryons and the form factors 
at interaction vertices but also the change of the internal (quark) structure of baryons in 
dense matter. The RHF calculations have performed in two ways: one with the coupling 



constants determined by SU(6) symmetry, 
on SU(3) (flavor) symmetry (see also Ref. 



;he other with the coupling constants based 



111]). Then, we have found that the baryon 



composition of the core matter in SU(3) symmetry is completely different from that in 
SU(6) symmetry. In SU(6) symmetry, all octet baryons appear in the density region below 
~ 1.2 fm~^, while, in the SU(3) calculation, only the S~ hyperon is produced. Furthermore, 
the medium modification of the internal baryon structure hardens the EoS for the core. 
Taking all those effects into account, we have obtained the maximum mass of a neutron star 
which is consistent with PSR J1614-2230. Therefore, it is very vital to consider the extension 
from SU(6) symmetry to SU(3) symmetry and the effect of the internal baryon-structure 
variation in nuclear matter. 

Now it is interesting to construct the EoS based on SU(3) symmetry in RMF approxima- 
tion, and to see how the symmetry extension affects the EoS, because the RMF calculation 
is practically much simpler than the RHF one, and many people have thus proposed many 
useful RMF models, some of which are accurately calibrated by various experimental data 
not only on infinite nuclear matter but also on finite nuclei. In this paper, we extend several 
popular RMF models, such as the GMl, GM3, NL3, TMl, FSUGold and lU-FSU models, 
and study the properties of nuclear matter and the mass-radius relations of neutron stars 
using the isoscalar, vector-meson couplings to the octet baryons in SU(3) symmetry. We 
then compare the results in SU(3) symmetry with those in the (usual) SU(6) calculations. 

In addition, we propose RMF models including the variation of baryon structure in a 
dense medium. In such models, we also use the coupling constants determined in SU(3) 
symmetry, and compare the results with those calculated in SU(6) symmetry. To take 
the variation of the in-medium baryon structure into account, we use the quark-meson 
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coupling (QMC) [12|, Il3[ and the chiral quark-meson coupling (CQMC) [IJ] models. It is 



well recognized that the constituent quark mass in a hadron is generally given by the quark 
condensate, {qq). The quark mass (or {qq)) in nuclear matter may then be reduced from 
the value in vacuum, because of the condensed scalar (a) field depending on the nuclear 
density, namely the Lorentz-scalar, attractive interaction in nuclear matter. The decrease of 
the quark mass leads to the variation of baryon internal structure at the quark level. Such 
an effect is considered self-consistently in the QMC model. 

The CQMC model is an extended version of the QMC model, in which the quark-quark 
hyperfine interaction caused by the one-gluon exchange is included. In addition, the pion- 
exchange interaction based on chiral symmetry is also considered. The hyperfine interaction 
plays an impotent role in the baryon spectra in matter [14,115[. The QMC and CQMC models 



have been successfully applied in studying the properties of hadrons in nuclear matter |16| . 



finite nuclei 
Ref. y.) 
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18] , hypernuclei [19|, l20| and neutron stars 
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2l|. (For a review, see 



Using those models, we calculate the particle fractions, the meson fields and the EoS 
inside the core. Furthermore, we estimate the maximum mass and radius of a neutron star 
by solving the the Tolman-Oppenheimer-Volkoff (TOV) equation 23|, |2J]. In the present 



calculations, we also study the role of the strange mesons {a* and 0) in the EoS. In SU(3) 
symmetry, we then find that the models except for GM3, FSUGold and lU-FSU can explain 
the mass of PSR J1614-2230. In the GM3, FSUGold and lU-FSU models, although the 
maximum mass cannot reach 1.97±0.04Mq, the calculated mass is not far from that value. 
Therefore, the extension from SU(6) to SU(3) symmetry is very vital for sustaining a heavy 
neutron star. In addition, the strange vector-meson (0) and the medium variation of baryon 
structure also help prevent the collapse of a neutron star. 

In RMF models, various types of nonlinear potentials with respect to the meson fields 
are usually involved, and they are very significant to reproduce the saturation condition 
for symmetric nuclear matter and the properties of finite nuclei. Among them, especially 
the Cao;^ term hardens the EoS at high density, and thus enhances a neutron-star mass. 
Furthermore, the nonlinear isoscalar-isovector coupling, A^pu'^p'^, which is involved only in 
the FSUGold and lU-FSU models, plays a unique role in the particle fractions in the core. 
If the cr-S and o"*-E coupling constants are determined so as to fit the (repulsive) mean-field 
potential for the E in nuclear matter, the S hyperon usually tends to be excluded in the 
core of a neutron star. However, in the FSUGold and lU-FSU models, the S~ as well as the 
A and S~ can emerge with a considerable fraction even at rather low density, which may be 
caused by the A^pu'^p^ interaction. 

This paper is organized as follows. In Section [TTl a brief review for RMF models based 



on Quantum Hadrodynamics (QHD) 25| is presented. The usual RMF models, the QMC 
and CQMC models are then unified through the scalar polarizability. In Section IIIIl the 
SU(3) extension in the coupling constants of the isoscalar, vector mesons is explained. The 
parameters in various models are determined in Section HVl Numerical results and discussions 
are addressed in Section |Vl Finally, we give a summary in Section |Vll 



II. RELATIVISTIC MEAN-FIELD MODELS 

For describing the properties of the core of a neutron star, we extend the usual Lagrangian 
density in RMF approximation to include not only the a, u and p mesons but also the 
strange mesons, namely the isoscalar, Lorentz scalar (a*) and vector (0) mesons. The a* 
and mesons are predominantly comprised of ss quarks. Because the charge neutrality and 
/3 equilibrium conditions are imposed in the core, the leptons must be introduced as well. 
The Lagrangian density is thus chosen to be 



C = ^^B ilf^d^ - Mb (a, a*) - Qu^bIi,^^ - 9^Blf,(p^ - OpBlf^P^' ■ Ib 



B 



Ipl 



where 






W^j. = d^Uu - dyU^ , Pf,u = df,(t)y - dycf)^ , R^y = 8 ^pu ^ d^Pf, , (2) 



with ipBie.) the baryon (lepton) field, Ib the isospin matrix for baryon B and mi the lepton 
mass. The sum B runs over the octet baryons, A^ (proton and neutron). A, S+'°'~ and H°'^, 
and the sum d, is for the leptons, e~ and p~ . The cu-, 0- and p-B coupling constants are 
respectively denoted by (7^^, g^j^B and g^B- In Eq- dlD, Unl is a nonlinear potential, which is 
explained below. 

When the baryons are treated as point-like objects (as in QHD), the effective baryon 
mass, M^, in matter is simply expressed as 

M% (a, a*) = Mb - g^Bcr - ga*B<y* , (3) 

where Mb is the mass in vacuum, and g^^B and ga*B are the a- and a*-B coupling constants, 
respectively. We hereafter call the model in which the baryons are structureless the QHD- 
type model. 

In contrast, in the QMC and CQMC models, the coupling constants, g^^B and ga*B-i 
depend on the a and a* fields, which reflects the variation of internal baryon structure in 



matter [13|, 



22]. Such dependences are caused by the attractive interactions due to the a 



TABLE I. Values of as, bs, a!g and b'^ for the octet baryons in the QMC or CQMC model. We 
assume that the scalar, strange {a*) meson does not couple to the nucleon. 



B 



as (fm) 



QMC 



CQMC 



a's (fm) 



b' 



B 



as (fm) 



a's (fm) 



b'r 



B 



N 


0.179 


1.00 


— 


— 


0.118 


1.04 


— 


— 


A 


0.172 


1.00 


0.220 


1.00 


0.122 


1.09 


0.290 


1.00 


S 


0.177 


1.00 


0.223 


1.00 


0.184 


1.02 


0.277 


1.15 


;=■ 


0.166 


1.00 


0.215 


1.00 


0.181 


1.15 


0.292 


1.04 



and a* exchanges. Thus, the in-medium baryon mass can be written as 21 1 



M* ((7, a*) = Mb - gaB{o)o - g.*B{(r*)a* 



with the following, simple parameterizations 

9aB{cr) = QaBhB I 



y,u,u,y,u 

O-B / 



9a*B{<y 



9a*Bh'B 



1 - ^ {g^,^a* 



(4) 



(5) 
(6) 



where QaN and gcj*K are respectively the a-N and a*-K coupling constants at zero density. 
Here, we introduce four parameters, a^, 6^, i'b and 6^, for describing the mass, and their 
values are tabulated in Table [B The effect of the variation of baryon structure at the quark 
level can be described with the parameters a^ and a'^. In addition, in the CQMC model, 
the extra parameters , hn and h'p, are necessary to express the effect of hyperfine interaction 



between two quarks [ij. 
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20|. If we set a^ = and 6^ = 1, gasi'^) becomes identical to 



the a-B coupling constant in QHD. This is also true of the coupling (^^.^(cr*). 

In the QHD-type model, we add the following nonlinear (NL) potential to the Lagrangian 
density 



UNL{cr,u>',pn 



1 



Kp {^puj^) {pp ■ pn 



(7) 



so as to reproduce the measured properties of nuclear matter and finite nuclei, for example, 
the incompressibility of nuclear matter, K^, the density dependence of symmetry energy, 04, 
etc. Here, the potential involves four coupling constants, g2, ga, C3 and A^p{= Ayg'^j^g'^^). 



In RMF approximation, the meson fields are replaced by the constant mean-field values: 
a, u, a*, (f) and p (the p*^ field). The equations of motion for the meson fields in uniform 
matter are thus given by 



mla + g2(T^ + ^3^^ = J^ 9aBCBi(r)p% , (8) 

B 

ml,r = Y,9.'BC's{r)p%, (9) 

B 

{ml + 2K^pp^) Q + c-iOJ^ = Y^ g^BpB , (10) 

B 

^l^ = Y9<l>BpB , (11) 

B 

[ml + 2KpU^) p = J29pB{lBhpB , (12) 



where the scalar density, p%, and the baryon density, ps, read 



,-.^/"-,,,=^£i(^Li:) . (13) 

1 r'^- .. ,2 4, 



with /ci?^ being the Fermi momentum for baryon B. 

In Eqs. ([S]) and (jH]), C^ and C^ are respectively the scalar polarizabilities (or the scalar- 
density ratios) at the a-B and a*-B interactions. Here, the scalar polarizabilities is defined 
by the ratio of the scalar density of a confined quark field at finite density to that in vacuum. 
In the QMC or CQMC model, they can be expressed by the following parameterizations |9|, 
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m 



Cb(o-) = 6b [1 -aB(fi'aAfO-)] , (15) 

C'^{a*) = b'^[l-a'Ag^,j,a*)] , (16) 

where the parameters a^, Bb, cl'b and b'^ take the same values as in Eqs. (^ and (^ (see 
also Table [T]). In contrast, they become unity in the QHD-type model (recall a^ = a^ = 
and Lb = h'^ = 1). 



The total energy density, e, and pressure, P, in the core then read 



e = ^^ I ^^ dke[e + Mf{a, a' 



B -^0 

2—2 — 4 2 T2 2 —2 o a — 2 —2 

+ J2^r'dke[e + m^Y^' , (17) 

^ = 4T^f-), (18) 

where the total baryon density, hb, is given by a sum of each baryon density 

Ub = '^pB ■ (19) 

B 

In the QHD-type model, the pressure can also be expressed as 



^ 2-2 -'- -3 -'- -4 -'- 2-*2 

I 2-2,-'- -4,-'- 2 72,-'- 2-2 I A -2-2 

+ 2"^^^ + 4^3^^ + 2"^0<^ + 2"^pP + ^'^p^ P 

fcF^ ^4 






III. SU(3) SYMMETRY IN THE ISOSCALAR, VECTOR-MESON COUPLINGS 

To study the EoS and the properties of neutron stars, it is very interesting to extend 
SU(6) spin-flavor symmetry based on the quark model to the more general SU(3) flavor 



symmetry 



11 



26| . Restricting our interest to three quark flavors (up, down and strange), 
SU(3) symmetry can be regarded as a symmetry group of strong interaction. To consider 
combinations of the meson-baryon couplings, it is extremely useful to choose the SU(3)- 
invariant interaction Lagrangian. Using the matrix representations for the baryon octet, B, 
and meson nonet (singlet state. Mi, and octet state, Mg), the interaction Lagrangian can 
be written as a sum of three terms, namely one coming from the coupling of the meson 
singlet to the baryon octet {S term) and the other two terms from the interaction of the 

8 



meson octet and the baryons - one being the antisymmetric (F) term and the other being 



the symmetric {D) term 111. 1271. |28|: 



Li^t = -9sV2[aTT{[B,Ms]B)+{l-a)TT{{B,Ms}B)]-gi^TT{BB)TT{M,) , (21) 
where gi and Qs are respectively the couphng constants for the meson singlet and octet states, 



and a (0 < a < 1) is known as the F/{F + D) ratio. For details, see the references |27|, |28 |. 

We here focus on the isoscalar, vector-meson {u and (p) couplings to the octet baryons, 

because, as usual, the other coupling constants can be determined so as to reproduce the 



observed properties of nuclear matter and hypernuclei (as discussed in section 



CYlli. 



The 



physical u and (p mesons are described in terms of the pure singlet, |1), and octet, |8), states 
as 

a; = cos6'^, |1) + sin6'^„|8) , = - sin6'^ |1) + cos6'^, |8) , (22) 

with By being the mixing angle. 

In SU(3) symmetry, all possible combinations of the couplings are then determined by 
four parameters: the singlet and octet coupling constants, gi and gs, the F/{F + D) ratio for 
the vector mesons, ay, and the mixing angle, 6y. If we require the universality assumption 
for the (electric) F/{F + D) ratio, we find ay = 1 27|, |29|]. In the limit of the ideal mixing, 
the mixing angle is given by 

ei"^^"^ = tan-i ( ^ j ~ 35.26° . (23) 

Furthermore, if the coupling ratio, z, is chosen to be 

2 = ^ = ^ ~ 0.4082 , (24) 

91 Vq 



we can obtain the usual SU(6) relations: 



-^QujN = ^9ujA = -j^QwY. = 9lj3 , (25) 

2^2 , , 

25'0A = 25-02 = 9<i>E = -i^9iuN , 9<t>N = . (26) 



In the present calculation, we refer to the Nijmegen extended-soft-core (ESC) model [27 1 



to fix the mixing angle and z. At present, this model may be the most complete model for 



^ When SU(3) symmetry is applied to the isovector, vector mesons, the Fock term is, in fact, necessary to 



reproduce the observed symmetry energy 



io|. 



the baryon-baryon interaction. It can well describe not only the A^-A^ but also the Y-N and 
Y-Y interactions in terms of the meson exchanges based on SU(3) symmetry. This model 
has then suggested the values of Oy and z as 

e^ = 37.50° , z = 0.1949 . (27) 

We notice that the mixing angle is very close to the ideal value, while the value of z is much 
smaller than that in SU(6) symmetry. It may be expected that a small value of z helps 



enhance the coupling constants ll| . We can find the relations of the coupling constants in 
SU(3) symmetry as 

1 1-V^2;tan6'„ 

gwA - Qi^T. - — -7=— —guN , Qur. - — -7=— —guN , [^O) 

l + VJ-Ztant/^ l + ^^/6ZtaJLlUv 

VSz - tan 9y , , 

^*^ = 1 I A7 , — ^^-^ ' (29) 
1 + v3-2tan6't, 

— tan6'^ ^/Sz + tanOy , , 

94>A = a<t>T. = — — 1=— — —gujN , 5'ois = -— — 1=— — —gLuN ■ (30) 

1 + V 3-2 tan Uy 1 + y/6z tan Uy 

Therefore, once the value of gi_jN is given, the other coupling constants, (^^y and g^ps, are 
determined by Eqs. (|28D-(l30l). 



IV. MODELS 



We examine two types of RMF models. One is based on the QMC and CQMC models 



10 



a, 



2l|, in which the variation of internal baryon structure in matter is taken into account. 
In these models, it is not necessary to consider any NL potential for describing the properties 
of nuclear matter around the saturation density, n%. The other is the QHD-type models 



with t 
GM3 



le NL potential given in Eq. ([7]). In fact, we adopt the parameterizations of the GMl, 







NL3 [30|, TMl J31|, FSUGold |32| and lU-FSU [33] models. Some of those models 
are very popular, because they are accurately calibrated by using various experimental data 
on infinite nuclear matter and finite nuclei. 

A. SU(6) symmetry 

In the case of QHD-type, the coupling constants, g^N, g^jN and gpN, are determined so 
as to reproduce the binding energy per nucleon, wq, and symmetry energy, 04, at n%. The 

10 



parameters, g2, gs, C3 and A^p, in Eq. ([7]) are chosen to be the values given in the original 
papers. For the vector-meson couplings to hyperons, we use the SU(6) relations given in 
Eqs. (^5^ and ( 12B|1 . and the following coupling relations 

1 . , 

9pN = ^QpT. = gps , fi-pA = . (31) 

Furthermore, assuming that the a* meson does not couple to a nucleon {ga*N = 0), the 
couplings of a-Y and a*-Y may be determined as follows. In RMF approximation, the 
potential for hyperon Y in symmetric nuclear matter, Uy , may be calculated as 

U^^^ = -QaYO- + QojYi^ ■ (32) 

Thus, we can determine the coupling constants, g^y, if we take the following values suggested 
from the experimental data of hypernuclei: f/^ = —28 MeV, f/^ = +30 MeV and 

f/4^) = -18 MeV 



34- 



In addition, if we consider the Nagara event |37j , which may suggest that the depth of the 
potential between two As is about —5 MeV, we may be able to fix the coupling constant, 
5'o-*A) by assuming that f/^ ^ —5 MeV, where Uj^ is the potential for A in A-hyperon 
matter: 

f/f ) = -<7.Aa(^) - g^,^a<^^ + g^ACO^^^ + g^^^^^^ . (33) 

Here, the superscript (Y) stands for a quantity in y-hyperon matter. Furthermore, we 
assume the relation, g^*^; = Qa'A, which is presented by SU(6) symmetry, and determine the 
coupling constant, ga*E, using the relation U~^ ~ 2f/| 
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In the QMC and CQMC models, the NL interaction is not necessary and the coupling 
constants can be determined by the same way as in the QHD-type model. We, however, 
notice that, in Eqs. (132|1 and (133|1 . the coupling constants for the scalar mesons should be 
replaced by the field-dependent ones (see Eqs. (|5]) and IQ). 

In Tables [TTl HIH and llVt we list the coupling constants in SU(6) symmetry and the 
properties of symmetric nuclear matter at n^. 

B. SU(3) symmetry 

As discussed in section UTTj because the pure singlet- and octet-states are mixed in SU(3) 
symmetry, the meson as well as the a and u mesons contributes to the nuclear saturation 

11 



properties. Thus, we have to readjust the couphng constants to satisfy the saturation condi- 
tion, namely the binding energy per nucleon, wq, at n^. We suppose that g^'N = and the 
couphng constant, Qo-n, takes the same value as in SU(6) symmetry. The coupling constant, 
QpN, is fixed so as to reproduce the value of symmetry energy given in the original paper. 

Firstly, we consider the QHD-type models. Assuming that the couplings, g2 and g^, in 
the NL potential takes the values given in the original paper, we can determine not only 
g^N and C3 but also g^jy so as to reproduce the same saturation condition as in the original 
paper. We here notice that, because g^pN is related to g^jN through Eq. ( l29|) . g^N is not free. 
Even in the case where the quartic term of the u field is not involved, namely C3 = 0, it is 
possible to reproduce the same saturation condition, because the 0-meson contributions to 
the energy density and pressure are quadratic (see Eqs. flT7|) and fl20|) ). and they have the 
same forms as in the w-meson contributions. For the vector-meson couplings to hyperons, 
we use the SU(3) relations given in Eqs. fl28|) and fl30l) . and the relations for the p meson, 
Eq. (|3Tl) . For the couplings of a-Y and a*-Y, we may be able to use the same procedure as 
in SU(6) symmetry. 

Next, in the QMC and CQMC models, we can also reproduce the same properties of 
nuclear matter as in SU(6) symmetry by only readjusting the coupling constant, g^jy (thus, 
g^N is also varied through Eq. (!29|) ). The other coupling constants may be determined by 
the same ways as in SU(6) symmetry. 

The results in SU(3) symmetry are also given in Tables HTl IIIII and IIVI 

V. NUMERICAL RESULTS AND DISCUSSIONS 
A. Properties of symmetric nuclear matter 

As seen in Tables [TTl Hill and IIVI the properties of symmetric nuclear matter are well 
reproduced in all the models. 

For the QHD-type models, the results calculated by GMl, GM3 and NL3 are presented 
in Table UTTt and those by TMl, FSUGold and lU-FSU are in Table ITVl In the former group, 
the NL potential involves the self-interaction terms of the a meson, while, in the latter group, 
in addition, the NL terms of the vector {u and p) mesons are taken into account. In the 
GMl, GM3 and NL3 models, the same saturation properties can be achieved in both SU(6) 

12 



TABLE II. Coupling constants and properties of symmetric nuclear matter in the QMC and CQMC 
models. We assume that ga*N = and go-*A = 5(t*e- The hadron masses are taken as follows: 
Mn = 939 MeV, Ma = 1116 MeV, Ms = 1193 MeV, Ms = 1318 MeV, m„ = 550 MeV, m^ = 783 
MeV, rup = 770 MeV, rria* = 975 MeV and rriff, = 1020 MeV. The saturation condition for 
symmetric nuclear matter is supposed to be wo = —15.7 MeV at n^ = 0.15 fm^^. The symmetry 
energy is taken to be 04 = 32.5 MeV at n^, and the slope parameter of the symmetry energy is 
denoted by l|j. 



vector sym. 



QMC 



SU(6) 



SU(3) 



SU(6) 



CQMC 



SU(3) 



9aN 

gujN 

9pN 

9<t>N 

9aK 

9aT. 

9aE 

9a* A 

9a*E 

K, (MeV) 
a4 (MeV) 
L (MeV) 



— 


coupling constants — 




8.28 


8.28 




8.50 


8.24 


7.98 




9.45 


4.38 


4.38 




4.29 


— 


-2.72 




— 


5.01 


6.09 




4.97 


2.45 


3.53 




3.24 


2.67 


4.83 




2.59 


1.09 


0.00^ 




2.62 


7.53 


5.19 




8.46 


— properties 


of symmetric 


nuclear matter - 


0.80 


0.80 




0.76 


280 


280 




302 


32.5 


32.5 




32.5 


88.7 


88.7 




90.7 



8.50 
9.14 
4.29 
-3.12 
6.11 
4.51 
4.84 
1.17 
5.80 

0.76 
302 
32.5 
90.7 



The symmetry energy, 04, is defined in terms of the 2nd derivative of the total energy with respect to 
the difference between proton and neutron densities, and the slope parameter, L, is then given by the 
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m- 



derivative of 04 with respect to the baryon density 
^ Because the tr-meson contribution in the QMC model already gives f/^ = — 8 MeV at nPg, the 

additional, attractive force due to the a* meson is not required. 
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TABLE III. Coupling constants and properties of symmetric nuclear matter in the GMl, GM3 and 
NL3 models. The relations, g^^'N = and ga*A = 9o-*i;, are assumed. For the NL3 model, we take 



rua = 508.194 MeV, m^ = 782.501 MeV and rup = 763.000 MeV [301. The other masses are the 



same as in Table HIl 



GMl GM3 NL3 

vector sym. SU(6) SU(3) SU(6) SU(3) SU(6) SU(3) 



coupling constants 



gaN 


9.57 


9.57 


8.78 


8.78 


10.217 


10.217 


52 (fm-1) 


12.28 


12.28 


27.88 


27.88 


10.431 


10.431 


53 


-8.98 


-8.98 


-14.40 


-14.40 


-28.885 


-28.885 


guiN 


10.61 


10.26 


8.71 


8.43 


12.868 


12.450 


gpN 


4.10 


4.10 


4.27 


4.27 


4.474 


4.474 


g<t>N 


— 


-3.50 


— 


-2.88 


— 


-4.250 


gak 


5.84 


7.25 


5.32 


6.51 


6.269 


7.853 


guT. 


3.87 


5.28 


2.85 


4.04 


4.709 


6.293 


ga-E 


3.06 


5.87 


2.83 


5.20 


3.242 


6.408 


ga*A 


3.73 


2.60 


2.03 


1.95 


5.374 


4.174 


ga*s 


9.67 


6.82 


7.89 


5.55 


11.765 


8.378 




— 


properties of symmetric 


nuclear matter — 


- 




n% (fm-3) 


0.153 


0.153 


0.153 


0.153 


0.148 


0.148 


Wo (MeV) 


-16.3 


-16.3 


-16.3 


-16.3 


-16.299 


-16.299 


M*^/Mn 


0.70 


0.70 


0.78 


0.78 


0.60 


0.60 


K^ (MeV) 


300 


300 


240 


240 


271.76 


271.76 


04 (MeV) 


32.5 


32.5 


32.5 


32.5 


37.4 


37.4 


L (MeV) 


93.9 


93.9 


89.7 


89.7 


118.0 


118.0 



and SU(3) cases, as explained in section [TV Bl (see also Eq. ( ITOj) and ( fTTT) ). In contrast, in 
the TMl, FSUGold and lU-FSU models, the symmetry energy, incompressibility and slope 
parameter, L, in SU(3) symmetry are slightly changed from the original values (given in 
SU(6) symmetry), because the terms of C3 and/or K^p in the NL potential, which has the 
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TABLE IV. Coupling constants and properties of symmetric nuclear matter in the TMl, FSUGold 
and lU-FSU models. The relations, grj*N = and qn-*A = 5(t*e, are assumed. For the TMl model, 
we take Mn = 938 MeV and m^ = 511.198 MeV [31], while, for the FSUGold and lU-FSU models. 



m 



33(. The other masses are 



m^ = 491.500 MeV, m^ = 782.500 MeV and mp = 763.000 MeV 
the same as in Table HH 

TMl FSUGold lU-FSU 

vector sym. SU(6) SU(3) SU(6) SU(3) SU(6) SU(3) 

— coupling constants — 



QaN 


10.029 


10.029 


10.592 


10.592 


9.971 


9.971 


92 (fm-1) 


7.233 


7.233 


4.277 


4.277 


8.493 


8.493 


53 


0.618 


0.618 


49.856 


49.856 


0.488 


0.488 


C3 


71.308 


81.601 


418.394 


522.820 


144.220 


171.586 


A ^ 


— 


— 


212.427 


168.100 


360.714 


248.010 


guiN 


12.614 


12.199 


14.302 


13.874 


13.032 


12.615 


9pN 


4.632 


4.640 


5.884 


5.395 


6.795 


5.821 


34>N 


— 


-4.164 


— 


-4.736 


— 


-4.306 


gak 


6.170 


7.733 


6.501 


8.295 


6.090 


7.680 


gaT. 


4.472 


6.035 


4.820 


6.615 


4.517 


6.107 


gas 


3.202 


6.328 


3.366 


6.953 


3.154 


6.334 


ga*A 


5.015 


3.691 


5.994 


4.458 


5.476 


4.204 


ga*E 


11.516 


8.100 
- properties 


13.071 
of symmetric 


9.147 
Quclear matter — 


11.915 


8.437 


n% (fm-3) 


0.145 


0.145 


0.148 


0.148 


0.155 


0.155 


Wo (MeV) 


-16.3 


-16.3 


-16.30 


-16.30 


-16.40 


-16.40 


M*^/Mb 


0.634 


0.634 


0.61 


0.61 


0.61 


0.61 


K^ (MeV) 


281 


284 


230 


252 


231.2 


237.7 


a4 (MeV) 


36.9 


36.9 


32.59 


32.59 


31.30 


31.30 


L (MeV) 


110.9 


110.8 


60.3 


66.6 


47.2 


54.6 



The coupling constant, K^p, also varies, because it is defined by A^^p — ^vd^NdtN^ where A„ takes the 
same value in both SU(3) and SU(6) cases. 
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FIG. 1. Particle fractions, Yi, in the QMC and CQMC models (left: 



1.0 1.2 



CQMC). 



quartic dependence of the nuclear density in the energy density and pressure of matter, also 
take part in reproducing the saturation condition. 

Furthermore, we notice the following two points. Firstly, in the extension of SU(6) to 
SU(3) symmetry, the coupling, g^j^, becomes smaller, because the (total) repulsive force is 
attributed not only to the u but also to the 0, which is caused by the mixing in Eq. (l22|) . 
We note that the coupling constant, g(f,N, is negative, because the mean-field value of the 
meson has a negative sign (see Figs. Hl-E]). 

Secondly, the coupling, g^'Y, (or the a* field itself) in SU(3) symmetry is suppressed 
in all the models, comparing with that in the SU(6) case. In contrast, the a-Y couplings 
in SU(3) symmetry are more enhanced than in SU(6) symmetry. This enhancement may 
counterbalance the additional, repulsive force due to the meson in the Y-N interaction, 
because the (total) repulsive force in the SU(3) case is stronger than in the SU(6) case. 

B. Neutron Stars 



In the core of a neutron star, the charge neutrality and /3 equilibrium under weak processes 
are imposed in solving the TOV equation 23|, |2J| . To obtain the realistic relation between 
the mass and radius of a neutron star, for the EoS at very low nuclear densities (< 0.068 
fm~^), we use the models given by Baym, Bethe, Pethick and Sutherland 40|, |4l|. In fact, 
the radius is relatively sensitive to the EoS at low density. 
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FIG. 2. Particle fractions, Yi, in the GMl, GM3, NL3 and TMl models (upper left: GMl, upper 
right: GM3, lower left: NL3, lower right: TMl). 



In Figs. [I]-[3l we show the particle fractions in the core of a neutron star. In each model, 
we calculate three cases: (1) only the non-strange mesons {a, u and p) are included in SU(6) 
symmetry; (2) all the mesons including the a* and are considered in SU(6) symmetry; (3) 
all the mesons are included in SU(3) symmetry. We note that, in some panels in the figures, 
the calculation stops at a certain density because the effective nucleon mass becomes zero 
beyond that density. 

As seen in the figures, from case (1) to (3) in order, the hyperons are created at higher 
densities. For example, the threshold densities of the A and S^ productions in SU(3) 
symmetry are higher than those in SU(6) symmetry, which makes the fractions of hyperons 
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FIG. 3. Particle fractions, Yi, in the FSUGold and lU-FSU models (left: FSUGold, right: lU-FSU). 




FIG. 4. Meson fields in the QMC and CQMC models (left: QMC, right: CQMC). 

small at high densities and thus increases the neutron fraction. 

In the models except for FSUGold and lU-FSU, the A and H"'^ hyperons are created, 
but the S does not appear, because the S-hyperon potential in nuclear matter, f/^ ; is 



chosen to be repulsive (see section lIVAp l. However, in the FSUGold and lU-FSU models, 



because the rather strong, w-p (nonlinear) repulsive interaction, K^pCj p % is included (see 
Eq-dZD), the H°'~ fields are very suppressed at high densities (especially, in the SU(3) case), 
and the S~ alternatively appears beyond ub — 0.4 — 0.6 fm~^. Furthermore, the order of 
the threshold densities for the S~ and S~ is reversed in the SU(6) and SU(3) cases. This is 

m 



^ We note that, if the Fock term is included 



42| . the hyperons except the ^ disappear. 
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FIG. 5. Meson fields in the GMl, GM3, NL3 and TMl models (upper left: GMl, upper right: 
GM3, lower left: NL3, lower right: TMl). 



a very remarkable phenomenon, and the isoscalar-isovector nonlinear interaction, K^pU p\ 
plays a unique role in the particle fractions. 

The meson fields are presented in Figs. |l]-|6l As it should be, in SU(6) symmetry, the 
strange-meson fields appear in the density region where the hyperons are generated. On the 
other hand, in SU(3) symmetry, the meson contributes to the baryon interactions even 
at low densities because of the mixing effect. However, the a* meson emerges above the 
density at which the first hyperon (usually the A) is created. This is because we assume 
that ga*N = 0. In the FSUGold and lU-FSU models, the fields of co and p mesons (especially 
the p) are very suppressed because of the isoscalar-isovector, nonlinear interaction. 
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FIG. 6. Meson fields in the FSUGold and lU-FSU models (left: FSUGold, right: lU-FSU). 
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FIG. 7. Equations of state in the QMC and CQMC models. 

In Figs. [7]- El we show the EoS in each model. Furthermore, in Figs. [10]- [121 we present 
the mass-radius relation of a neutron star calculated by the TOV equation. The detail of 
the neutron-star properties is also shown in Table [Vl 

As expected, because the isoscalar, vector-meson couplings to the octet baryons are en- 
hanced in SU(3) symmetry, the extension from SU(6) to SU(3) symmetry hardens the EoS 
very much. In each model, the hardest EoS is given by the case (3), while the softest one is 
obtained in the case (1). This tendency can be related to the fact that, as seen in Figs. [Tj- 
[31 the densities at which the hyperons appear in the case (3) are rather higher than those in 
the case (1). In general, the strange mesons, especially the meson, also play an important 
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FIG. 8. Equations of state in the GMl, GM3, NL3 and TMl models (left: GMl and GM3, right: 
NL3 and TMl). 
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FIG. 9. Equations of state in the FSUGold and lU-FSU models. 



role in supporting a heavy neutron star. In the mass-radius relations presented in Figs. [TO] - 
T2I we can again see that, in each model, the maximum neutron-star mass in the case (3) is 
heaviest, while the lightest one is given in the case (1). We note that, in Fig. [TTl the curve 
(red solid) for the NL3 model in SU(3) symmetry cannot reach the maximum point because 
the nucleon mass becomes negative at the density before the maximum point. 

We summarize the following, several comments on the mass-radius relations shown in 
Figs, [in]- [121 In the QMC and CQMC models, the maximum neutron-star masses calculated 
in SU(3) symmetry are consistent with the pulsar PSR J1614-2230. In particular, the mass 
in the CQMC model clearly exceeds the mass of PSR J1614-2230. Because the difference 
between the QMC and CQMC models is originated by the hyperfine interaction between 
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FIG. 10. Mass-radius relations in the QMC and CQMC models. 
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FIG. 11. Mass-radius relations in the GMl, GM3, NL3 and TMl models (left: GMl and GM3, 
right: NL3 and TMl). In the NL3 model, the (red) solid curve for SU(3) symmetry does not yet 
reach the maximum point because the nucleon mass in matter becomes negative at the endpoint. 

two quarks inside a baryon, the large difference between the two maximum masses is mainly 
generated by this microscopic interaction. It is noticeable that the quark-quark hyperfine 
interaction is very vital to obtain the correct mass spectra of octet baryons in a nuclear 
medium [lJ,|2o|. 

In the GMl model, the maximum neutron-star mass in SU(3) symmetry is much larger 
than the observed mass of PSR J1614-2230. In contrast, the maximum mass in the GM3 
model is clearly under the observed value (see the left panel in Fig. [TTll . The difference 
between the two models is just in the values of the nuclear incompressibility and the slope 
parameter, namely K^ = 300 (240) MeV and L = 93.9 (89.7) MeV for the GM1(3) model 
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FIG. 12. Mass-radius relations in the FSUGold and lU-FSU models. 



(see Table [ITTD. 

The NL3 model is a unique model, in which the mass of PSR J1614-2230 can be explained 
even in SU(6) symmetry (see the right panel in Fig. [TT]) . This model may be characterized 
by the large values of symmetry energy (04 = 37.4 MeV) and slope parameter (L = 118 
MeV) (see Table [ml). 

The rather large values of 04 and L are also used in the TMl model, where only the SU(3) 
result can, however, reach the mass of PSR J1614-2230. Furthermore, the difference between 
the maximum masses in SU(6) and SU(3) symmetries is very large in the TMl model (see 
also Table |V|) . This fact may be caused by the repulsive force due to the nonlinear csu^ term 
in Eq.([7]). Note that, to reproduce the same saturation condition as in SU(6) symmetry, the 
strength of C3 in SU(3) symmetry is larger than that in SU(6) symmetry (see Table HVT) . 

Unfortunately, the mass of PSR J1614-2230 cannot be explained by the FSUGold and 
lU-FSU models. However, in both models the maximum mass in SU(3) symmetry becomes 
1.8 — 1.9Mq, which is not far from the observed mass. In these models, the maximum mass 
in SU(6) symmetry is again very different from the value in SU(3) symmetry (see Table IVT) . 
Furthermore, although the curves for M/Mq in the SU(6) and SU(3) cases normally coincide 
with each other in the low mass region (see Figs. [TU]- [TT]) . the two curves clearly stay away 
from each other even at M/Mq = 0.8 in the FSUGold and lU-FSU models (see Fig. [T^. 
These facts may again be caused by the very large difference between the values of C3 in 
SU(6) and SU(3) symmetries (see Table HVj) . 
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TABLE V. Properties of a neutron star in SU(6) or SU(3) symmetry. We list the neutron-star 
radius, i?max (in km), the ratio of the neutron-star mass to the solar mass, M^^-^/Mq^ and the 
central density, Uc (in fm~^) at the maximum-mass point. In these calculations, we consider all 
the mesons {a^ uj, p, a* and </>). 

SU(6) SU(3) 

^max Mmax/M0 Tic i?max Mmax/Mg Tic 



QMC 


12.5 


1.72 


0.85 


11.8 


1.93 


0.96 


CQMC 


12.6 


1.84 


0.84 


12.1 


2.08 


0.90 


GMl 


12.7 


1.86 


0.82 


12.2 


2.14 


0.87 


GM3 


12.1 


1.63 


0.93 


11.4 


1.85 


1.05 


NL32 


13.1 


2.07 


0.78 


— 


— 


— 


TMl 


13.1 


1.72 


0.77 


12.5 


2.03 


0.86 


FSUGold 


11.4 


1.39 


1.03 


11.2 


1.79 


1.08 


lU-FSU 


11.3 


1.55 


1.03 


11.3 


1.88 


1.02 



'^ In the NL3 model with SU(3) symmetry, the nucleon mass becomes negative before the neutron-star 
mass reaches the maximum point. Therefore, the maximum mass is not given in SU(3) symmetry. 

VI. SUMMARY 

We have calculated the particle fractions, the meson fields and the EoS in the core of a 
neutron star, using the popular RMF models (such as GMl, GM3, NL3, TMl, FSUGold 
and lU-FSU) as well as the QMC and CQMC models. It is noticeable that some of the RMF 
models are accurately parameterized to compute the properties of infinite nuclear matter 
and finite nuclei. On the other hand, because, in the QMC and CQMC models, the quark 
degrees of freedom in a baryon are taken into account, they allow us to consider the variation 
of the quark structure of baryon in dense mater. In particular, the CQMC model involves 
the quark-quark hyperfine interaction, and thus it can correctly describe the octet baryon 
spectra in matter as well as in vacuum. 

In the present calculations, we have examined the extension from SU(6) spin-flavor sym- 
metry based on the quark model to SU(3) flavor symmetry in determining the isoscalar, 
vector-meson coupUngs to the octet baryons. We have also studied how the strange mesons 
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{a* and 0) contribute to the internal structure of a neutron star. 

In SU(3) symmetry, we have found that the models except GM3, FSUGold and lU-FSU 
can explain the mass of PSR J1614-2230. In the GM3, FSUGold and lU-FSU models, 
although the maximum mass cannot reach 1.97 ± O.O4M0, the calculated mass is not far 
from that value. Therefore, the extension from SU(6) to SU(3) symmetry and the strange 
vector meson, 0, are very significant in sustaining a heavy neutron star. In addition, the 
variation of baryon structure in matter also helps prevent the collapse of a neutron star. 

In RMF models, the NL potential, U^l, is indispensable for reproducing the satura- 
tion condition for symmetric nuclear matter and the properties of finite nuclei. In the 
present calculations, it involves not only the usual, nonlinear a terms but also the csw^ term 
and the isoscalar-isovector A^pW^p^ term. Among them, in particular, the nonlinear c^oo^ 
term hardens the EoS, and thus enhances a neutron-star mass. Furthermore, the nonlinear 
isoscalar-isovector coupling plays a unique role in the particle fractions inside a neutron 
star. Because the a-T, and (t*-S coupling constants are usually determined so as to fit the 
(repulsive) mean-field potential for the S in nuclear matter, it becomes difficult to create 
the S hyperon in the core of a neutron star. However, in the FSUGold and lU-FSU mod- 
els, instead of the S°'~, the S~ can emerge with a considerable fraction even at rather low 



density. As the power counting [43|, |4J] suggests, there may be many possible, NL couplings 
and many-body forces containing various meson fields, which may contribute to the EoS. It 
is thus interesting to study how such interactions contribute to the properties of a neutron 
star. We note that, at the MF level, the NL potential may be regarded as many-body inter- 
actions among baryons because the meson fields are just auxiliary fields and thus they can 
be replaced with bilinear forms of baryon fields. 

In RMF models, the parameterization is usually performed using experimental data mea- 
sured around n^. However, because the region of density which is important in the EoS for 
a neutron star may be the region of 0.8 — LI fm^'^ (> Qn%), no one knows if such parame- 
terizations work correctly at such high densities. Therefore, although in this paper we have 
studied RMF models from various standpoints, it may be difficult to winnow the good mod- 
els out at the RMF level. As suggested in Dirac-Brueckner-Hartree-Fock calculations 45|, it 
may, at least, be imperative to include the density dependence of the parameters to obtain 
conclusive results on the EoS. 

In the present calculations, we have not considered the Fock (exchange) term. Although, 
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in naive QHD-I [25], the Fock contribution seems very small in symmetric nuclear matter, 
it plays a very important role even around the normal nuclear matter density as well as at 
high densities, if the p meson is included and the tensor interaction thus arises |9|, |lOJ. It is 
remarkable that, when the tensor interaction is taken into account, the A hyperon does not 
appear in the core even at high density |9|, [lO|, |42|. Furthermore, the tensor contribution is 



very important in reproducing the density dependence of symmetry energy, 04. 

At very high density, the quark and gluon degrees of freedom, rather than the hadron 



degrees of freedom, may take place in the core matter [46| . Because the degrees of freedom in 
quark-gluon matter are generally large, it is necessary to assume a rather strong correlation 
between quarks and gluons to support a massive neutron-star mass. It would be very 
interesting to investigate how the quark-gluon phase connects with the hadron one and how 
such degrees of freedom contribute to the EoS for a neutron star. 
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